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Using a Monte Carlo ray-tracing code in full general relativity, we calculate the transport of 
photons from a geodesic hot spot emitter through a corona of hot electrons surrounding a black 
hole. Each photon is followed until it is either captured by the black hole or is detected by a distant 
observer. The source is assumed to be a low-energy thermal emitter (T cm ~ 1 keV), isotropic in 
the rest frame of a massive geodesic test particle. The coronal scattering has two major observable 
effects: the Comptonization of the photon spectrum due to the high energy electrons, and the 
convolution of the time-dependent light curve as each photon is effectively scattered into a different 
time bin. Both of these effects are clearly present in the Rossi X-Ray Timing Explorer (RXTE) 
observations of high frequency quasi-periodic oscillations (QPOs) seen in black hole binaries. These 
QPOs tend to occur when the system is in the Steep Power Law spectral state and also show no 
evidence for significant power at higher harmonic frequencies, consistent with the smoothing out 
of the light curve by multiple random time delays. We present simulated photon spectra and light 
curves and compare with RXTE data, allowing us to infer the properties of the corona as well as 
the hot spot emitter. 



I. INTRODUCTION 

In recent years, observations of accreting black holes 
with the Rossi X-ray Timing Explorer (RXTE) have 
discovered a number of sources with high frequency 
quasi-periodic oscillations (HFQPOs) in their X-ray 
light curves. These sources are seen predominately in 
the Steep Power Law (SPL) spectral state, suggesting 
a significant level of inverse-Compton scattering of the 
emitted photons off of hot coronal electrons (for an 
excellent review of the observations, see 0). 

Motivated by these observations, we extend a sim- 
ple geodesic hot spot model [n|, [l2| to include Monte 
Carlo scattering of photons emitted isotropically in 
the hot spot's rest frame, which then propagate 
through a corona surrounding the black hole. The hot 
spot has a planar orbit near the inner-most stable cir- 
cular orbit (ISCO) and the corona is modeled with a 
self-similar density profile described by the advection 
dominated accretion flow (ADAF) model 

In Section [H] we describe the physics of classical 
electron scattering of unpolarized light, and show how 
energy is transferred from high-energy electrons to 
low-energy photons. Section iLUl shows how this scat- 
tering is treated in a relativistic context in the Kerr 
metric. Assuming a thermal emitter with T em = 1 
keV, we show in Section IIVI how the observed spec- 
trum is modified by the hot electrons. The shape of 
this modified thermal spectrum can in turn be used 
to infer the properties of the corona. 

Section Ivl shows the time-dependent effects of scat- 
tering on the X-ray light curve, particularly how the 
amplitude of modulation is damped for larger optical 
depths. Furthermore, we find this damping is more 
significant for photons experiencing multiple scatter- 
ing events, which tend to have higher energies. We 
conclude in Section IVII with a summary of the impli- 
cations these results have for the hot spot model for 



QPOs. 



II. ELECTRON SCATTERING 



Following Rybicki & Lightman [l(j , we use the cross 
section for Thomson scattering of unpolarized radia- 
tion off of nonrelativistic electrons: 
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where rg is the classical electron radius ro 



(1) 
2.82 x 



It is important to note that "nonrelativistic" is a ref- 
erence to the photon energy, not the electron energy. 
In the electron rest frame, we require hv <C m e c 2 in 
order for the above cross sections to be valid, in which 
case the scattering is nearly elastic or coherent. For 
higher energy photons, the scattering involves quan- 
tum effects and requires the "Klein-Nishina" cross sec- 
tion [j| . Since we are primarily interested in the scat- 
tering of photons from a relatively cool thermal ac- 
cretion disk {hv ~ 1 — 5 keV), the classical treatment 
should suffice. 

Even if the scattering is treated as coherent in the 
electron frame, in the lab frame energy can be (and 
often is) transferred from the electron to the photon. 
To see this boosting effect, consider a photon with 
initial energy Ei scattering off an electron with velocity 
(3 in the ^-direction in the "lab frame" K. In this 
frame, the angle between the incoming photon and 
electron velocity is 9. In the electron rest frame K', 
the photon is scattered at an angle 9' with respect to 
the x'-axis. The Doppler shift formula gives 

e- = Sij(l — (3 cos 9) 



£y7(l + /3cos 6'), 



(2) 



where 7 = l/yl — P 2 and £/ is the post-scattering 
energy in the lab frame. In the electron frame, we 
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assume elastic scattering with e[ = e'p which should 
be the case for the typical seed photons from a thermal 
emitter at T cm ~ 1 keV. 

Averaging over all angles 6 (isotropic electron dis- 
tribution) and 6' [weighted by eqn. we find that 
the typical scattering event boosts the photon energy 
by 



Si 



7 



(3) 



To determine 7, we consider a Maxwell-Boltzmann 
distribution function in electron momentum p = r ymv: 



f(p)d 3 p oc 4irp 2 exp — 



Vp 1 



kT P _ 



(4) 



III. RELATIVISTIC IMPLEMENTATION 

Unlike the approach taken in [llj, where the pho- 
tons were traced backwards in time from a distant ob- 
server to the emitting region, here it is conceptually 
easier to trace the photons forward in time from the 
emitter to the observer, then use Monte Carlo meth- 
ods to determine the distribution of scattered photons. 
For the Thomson cross section, the path of each pho- 
ton is energy- independent, so the photon's final ob- 
served energy can be thought of as a fiducial redshift 
E y> s /E cm that can be convolved with the spectrum in 
the local emitter frame to produce the total spectrum 
seen by the observer. 

To determine the initial momentum of each photon, 
we construct a tetrad centered on the emitter's rest 
frame, denoted by tilde indices fi. In the coordinate 
basis, the energy and angular momentum of a particle 
on a stable circular orbit around a Kerr black hole are 
given by 



-pt(em) = 



r 2 - 2Mr ± a^/Wr 
r(r 2 ~ 3Mr ± 2a v / Mr) 1 / 2 



and 



i (em) = ± 



VMr(r 2 =p 2«VMr + a 2 ) 
r(r 2 - 3Mr ± ZaVMr') 1 / 2 



(5) 



(6) 



From these we construct the 4- velocity via the inverse 
metric p tJ, (em) — g fiU p u (em), which gives e t ~. Then 
ef and are defined parallel to the coordinate basis 
vectors e r and eg, and is given by orthogonality. 

In this basis, the initial photon direction is picked 
randomly from an isotropic distribution, uniform in 
spherical coordinates cos 6 = [—1,1] and <\> = [0, 2ir). 
All photons are given the same initial energy in the 
emitter frame pi = —E Q , which is used as a reference 
energy for calculating the final redshift with respect 
to a stationary observer at infinity. Given p^, the 




FIG. 1: Planar photons emitted isotropically in the rest 
frame of a massive particle on a circular orbit at the ISCO. 
The photon paths are colored according to their red- or 
blue-shift in energy with respect to Eo measured in the 
emitter's frame. 



photon's geodesic trajectory is integrated using the 
Hamiltonian formulation described in |ll| . 

Figure ^ shows an "overhead view" of photon tra- 
jectories in the plane of the disk, emitted isotropically 
by a massive test particle on a circular orbit at the 
ISCO of a black hole with a/M = 0.5. The pho- 
tons are colored according to their energy-at-infinity 



E n 



-pt, either blue- or redshifted with re- 



spect to their energy in the emitter frame Eq>. The rel- 
ativistic beaming is done automatically by the Lorentz 
boost from the emitter to the coordinate frame, so 
the blue photons are clearly bunched more tightly to- 
gether, as required by the invariance of I u /v 3 . 

At each step along the photon's path, we determine 
the probability of electron scattering according to the 
differential optical depth dr cs — n cs pds. The density p 
is defined in the "Zero Angular Momentum Observer" 
(ZAMO) frame and the opacity k cs is given by 
the classical cross section quoted above in equation 
JTJ. The scattering is then treated classically in the 
electron's local frame, after which the new photon 4- 
momentum is transformed back to the coordinate ba- 
sis and the ray-tracing continues along the new tra- 
jectory. 



IV. EFFECT ON SPECTRA 

As we showed at the end of the Section [D] one ef- 
fect of the coronal scattering is generally a transfer of 
energy from the electron to the photon. One way to 
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quantify this energy transfer is through the Compton 
y parameter, defined as the average fractional energy 
change per scattering, times the number of scatterings 
through a finite medium. For nonrelativistic electrons, 
Rybicki & Lightman ^(| show that the average energy 
transfer per scattering event is 



1 o 5 r 



AkT„ 



m e c* 



(7) 



The mean number of scatterings for an optically 
thin medium is simply r cs , the total optical depth 
through the medium. For optically thick systems, the 
photons must take a random walk to escape, so the 
number of scatterings becomes t 2 s . Thus the Comp- 
ton y parameter for a finite medium of nonrelativistic 
electrons is 



AkT R 



m e c^ 



Max(r< 



'cs)- 



(8) 



For a low-energy soft photon source with multiple 
scattering events, the final spectrum due to inverse- 
Compton scattering can be estimated using the Kom- 
paneets equation. For kT em < hv < kT e , the resulting 
spectrum takes the power-law form 



with 



(9) 



(10) 



At energies above kT e , the electrons no longer effi- 
ciently transfer energy to the photons, so the spectrum 
shows a cutoff for hv > kT e : 



I y ~ v 3 exp(— hv/kT e ). 



(11) 



With the assumption of purely elastic scattering, we 
cannot actually reproduce this cutoff effect; all pho- 
tons are scattered equally, and the ratio £//£» is inde- 
pendent of energy. Thus equation Q would predict 
infinite energy boosts until hv m e c 2 . In reality, 
higher energy photons tend to lose energy in scatter- 
ing, due to the recoil of the electron. This effect is 
relatively easy to calculate from conservation of en- 
ergy and momentum in the electron rest frame: 



1 + -^s(l - cos6>) 



(12) 



To accurately include this effect, we would have to 
keep track of the real "physical" energy of each pho- 
ton, instead of the fiducial redshift method that we 
currently use to reconstruct the total spectrum after- 
wards. Ultimately, this is just a matter of computa- 
tional intensity and no real conceptual difficulty. To 
first-order, we can treat the thermal photon source as 
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FIG. 2: Simulated observed spectra of a thermal hot spot 
emitter with T cm = 1 keV, on a circular orbit at the ISCO 
of a black hole with a/M — 0.5. The thermal spectrum is 
modified by relativistic effects and Compton scattering off 
a hot corona with T e = (nsco/f)100 keV. 



a monochromatic emitter at Eq — 3fcT cm , which gives 
a reasonable approximation to the true solution. 

For the corona geometry, we use a self-similar distri- 
bution described by an ADAF model 8] , with density 
and temperature profiles that scale as 

-3/2 



p OC 

T oc 



(13) 
(14) 



outside of the ISCO. We have ignored the bulk velocity 
of the inwardly flowing gas, which will typically have 
fbuik <C fthcrm in the ADAF model. 

Without scattering, the time-averaged "numerical" 
spectrum could be described by the relativistic trans- 
fer function described in |ll| , defined over an infinites- 
imal band in radius i?; n ~ Rout — f cm . The inverse- 
Compton processes in the corona serve to further 
broaden this transfer function, as in j3| and . This 
transfer function is then normalized to the rest energy 
Eq and convolved with the actual emission spectrum 
(e.g. a thermal blackbody at kT em ) to give the simu- 
lated observed spectrum. 

Figure [21 shows a set of these simulated spectra 
from a hot spot emitter around a black hole with 
a/M — 0.5. The emission spectrum is thermal in 
the hot spot rest frame with T em — 1 keV. The coro- 
nal ADAF model has T e = (ri SCO /r)100 keV, and 
electron density n e ~ r ~ 3 ' 2 for a variety of opti- 
cal depths t cs . The spectra are plotted in units of 
(^/=photons/s/cm 2 /keV), as is the convention by many 
observers, but the actual magnitude of the y-axis is ar- 
bitrary, and would normally depend on the distance 
to the source. From the slope of the power-law and 
the location of the cut-off, the corona temperature 
and optical depth can be inferred from observations 

aula 
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V. EFFECT ON LIGHT CURVES 

The spectra in Figure [21 were created by integrating 
over the complete hot spot orbital period and over 
all observer inclination angles. However, during the 
Monte Carlo calculation, it is just as simple to bin all 
the photons according to their final values of 9, t bs) 
and energy —pt- The latitude bins are evenly spaced 
in cos so that a comparable number of photons land 
in each zone. The energy bins are spaced logarithmi- 
cally to include the high energy tail and also maintain 
high enough spectral resolution at lower energies. As- 
suming the hot spot is on a circular periodic orbit, 
the photons detected at any azimuthal position can 
be mapped into the appropriate bin in i b s , so that 
none are "wasted." 

An excellent way to see the effects of scattering 
on the hot spot light curves is by plotting time- 
dependent spectra of a monochromatic emitter, shown 
in Figure [3] for a range of optical depths and an in- 
clination angle of 75°. The logarithmic color scale 
shows (#photons/s/cm 2 /keV/period), normalized to 
the peak intensity in each panel. At "0" phase, when 
the emitter is on the far side of the black hole, the 
spectrum shows two distinct lines, one blueshifted in 
the forward direction of hot spot motion, and one 
redshifted in the backward direction. As the hot 
spot comes around towards the observer, the directly 
beamed blueshifted line dominates, and then when the 
phase is ~ 0.5 and the hot spot is on the near side of 
the black hole, a single line dominates. This is due 
to the gravitational demagnification of the secondary 
images formed by photons that have to complete a 
full circle around the black hole to reach the observer. 
While these features would most likely be unresolvable 
for black hole binaries, they may well be observable in 
X-ray flares from Sgr A* as well as other supermassive 
black holes (e.g. see pj). 

In the subsequent panels, the spectrum is modified 
by the scattering of the hot spot photons in the sur- 
rounding corona. As in Section IIVI the temperature 
and density profile of the corona is given by an ADAF 
model with T e (nsco) = 100 keV. The four panels 
of Figure show increasing values of T es = [0, 1, 2, 4]. 
The effects of scattering on the spectra are really quite 
profound. As we described qualitatively in |l2(, the 
electron corona is like a cloud of fog surrounding a 
lighthouse, spreading out the delta-function beam in 
time and energy. 

Unlike the simple model there, where each pho- 
ton was assigned some positive time delay, the full 
Monte Carlo scattering calculation shows that some 
photons actually arrive earlier in time by taking a 
"shortcut" to the observer instead of waiting for the 
hot spot to come around and move towards the ob- 
server. And of course, the photons are also spread 
out in energy due to the inverse-Compton effects. As 
the optical depth increases, the well-defined curves in 



Figure |3Ji is smeared out into a nearly constant blur 
when r cs = 4, with a broad spectral peak as in Figure 
13 Only a slight trace of the original coherent light 
curve remains, composed of roughly 1% of the emit- 
ted photons that do not scatter before reaching the 
observer or get captured by the event horizon. When 
r es > 1, multiple scattering become more common, 
so photon shortcuts become rarer, tending to spread 
the light curve preferentially to the right (delay in ob- 
server time), as seen in Figures |3J;,d. 

By integrating over broad energy bands such as 
those typically used in RXTE observations, we can 
increase our "signal" strength while sacrificing spec- 
tral resolution. For millisecond periods, there will still 
not be nearly enough photons to provide phase reso- 
lution, but these features may show up statistically 
in the power spectrum or bispectrum |(| . Figure 0] 
shows a set of integrated light curves for a variety of 
optical depths. The black hole and hot spot parame- 
ters are as in Figure |3 here assuming a thermal hot 
spot emitter with temperature T^ s = 1 keV. As the 
optical depth to electron scattering increases, the rms 
amplitude of each light curve decreases as the pho- 
tons get smoothed out in time. Similarly, due to the 
average time delay added to each photon by the in- 
creased path length, the relative location of each peak 
is shifted later in time. 

In all likelihood, the relative phase shifts would be 
nearly impossible to detect, regardless of the instru- 
ment sensitivity, since to do so would require measur- 
ing the light curve from a single coherent hot spot at 
two different optical depths. It is difficult to imagine 
a scenario where the coronal properties could change 
on such short time scales (yet it is possible that a fixed 
hot spot on the surface of an X-ray pulsar might actu- 
ally be used for this technique). However, the higher 
harmonic peaks of the different light curves may in 
fact be measurable with the next-generation X-ray 
timing mission, or under extremely favorable condi- 
tions, even with RXTE. 

While the absolute peak shifts for hot spot light 
curves at different optical depths would probably not 
be detectable, the relative shifts of simultaneous light 
curves in different energy bands may be observable, 
at least on a statistical level with a cross-correlation 
analysis. Since the average scattering event boosts 
photons to higher energy bands and also causes a net 
time delay due to the added geometric path, the light 
curves in higher energy bands should be delayed with 
respect to the lower energy light curves. A few of the 
typical energy bands used for RXTE observations are 
(2-6), (6 - 15), and (15 - 30) keV. To fully cover 
the peak emission from a thermal hot spot at 1 keV, 
we expand the lowest energy band in our calculations 
to cover (0.5 — 6) keV. The light curves in these three 
bands are plotted in Figure [5] for i = 75°. The low en- 
ergy band resembles the unscattered light curve plus a 
roughly flat background, while the higher energy light 
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FIG. 3: Time-dependent spectra of a monochromatic, isotropic hot spot emitter on an ISCO orbit with a/M = 0.5 and 
inclination angle i = 75°. The four panels show spectrograms for systems of increasing optical depth r cs = [0, 1,2,4]. 
The logarithmic color scale shows the number of photons in each time/energy bin, normalized to the peak value for each 
panel. 
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FIG. 4: Energy-integrated light curves for a hot spot with 
orbital parameters as in Figure [3] and a range of coronal 
optical depths. The emitted spectrum is assumed to be 
thermal with a hot spot temperature Th s = 1 keV, inte- 
grated over 0.5 — 30 keV in the observer's frame. 



curves show a much smaller modulation with a signif- 
icant phase shift (~ 0.3 periods) due to the additional 
photon path lengths. 

VI. IMPLICATIONS FOR QPO MODELS 

The original motivation for the application of scat- 
tering to the hot spot model was to answer a few im- 
portant questions raised by RXTE observations: 

• The distinct lack of power in higher harmonics 
at integer multiples of the peak frequencies. 

• The larger significance of high frequency QPO 
detections in the higher energy bands (6-30 keV) 
relative to the signal in the lower energy band 
(2-6 keV). 



• The trend for these HFQPOs to exist predom- 
inantly in the Steep Power Law (SPL) spectral 
state of the black hole. 

Beginning with the final point, it appears to be 
quite reasonable that the physical mechanism pro- 
ducing the power law region of the spectrum is the 
invcrse-Compton scattering of cool, thermal photons 
off of hot coronal electrons. The steep power law sug- 
gests a small-to-moderate value for the Compton y 
parameter, inferred from equation (|10l) . in the range 
0-5 < y < 10. From equation JSJ), this suggests either 
a small optical depth or a small electron temperature. 
To gain insight into which of these two options is more 
likely, we need to address the other two observational 
clues. 

In ^3], we nrs t proposed the scattering model as 
an explanation for harmonic damping. With the more 
careful treatment in this paper, we include not only 
the temporal, but also the spatial effects of electron 
scattering. The photons originally beamed toward the 
observer are now scattered in the opposite direction, 
while the photons emitted away from the observer can 
now be scattered back to him. This smoothes out the 
light curve in time more effectively than the localized 
convolution functions used in ^^l- At the same time, 
the scattering is not completely isotropic [see eqn. JQ], 
so some modulation remains. Thus, to maintain a 
significant modulation in the observed light curve, we 
require a relatively small optical depth, reducing the 
smoothing effects of the scattering. 

The fact that most HFQPOs appear more signif- 
icantly in higher energy bands also points towards 
Compton scattering off hot electrons. However, as 
the calculations above show (see Fig. [5}, with the basic 
thermal disk/hot spot model, the light curves actually 
have smaller amplitude fluctuations in the higher en- 
ergy bands, as these scattered photons get smoothed 
out more in time. Furthermore, while the higher har- 
monic modes are successfully damped in the scattering 
geometry, so is the fundamental peak. Thus, in order 
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FIG. 5: Hot spot light curves in a few different RXTE 
energy bands (we have expanded the lowest energy band 
down to 0.5 keV to include the thermal emission of a hot 
spot at Th s = 1 keV). The hot spot inclination is 75° and 
the coronal properties are as in Figure [3] The optical 
depth to scattering is r cs = 1.5 in (a) and 2.5 in (b). The 
higher energy light curves are made from photons that 
have experienced more scattering events, boosting their 
energy and delaying their arrival time. 



Based on these arguments alone, we find it unlikely 
that the HFQPOs are coming from a cool, thermal 
hot spot getting upscattered by a hot corona. From 
the photon continuum spectra of the SPL state, there 
appears to be a hot corona with T e ~ 100 keV, but 
as Figure shows, the lowest energy photons, which 
presumably do not scatter in the corona, have by far 
the greatest amplitude modulations. It is possible 
that the relative modulation would appear smaller due 
to the added flux from the rest of the cool, thermal 
disk, but much of this steady-state emission should 
also get scattered to higher energies, further damping 
the modulations in the (6 — 30) keV bands. 

The high luminosity of the SPL state (also called 
the Very High state) suggests that the thermal, slim 
disk geometry may not be appropriate here. Per- 
haps it is more likely that these cases correspond to 
an ADAF model, traditionally associated with very 
low or very high accretion rates Q . Since the ADAF 
model cannot radiate energy efficiently, the gas in the 
innermost regions will be much hotter than in the slim 
disk paradigm. Thus hot hot spots with Th s > 5 keV 
could be forming inside a small ADAF coronal region, 
providing seed photons that are already in the higher 
energy bands, and are only moderately upscattered by 
the surrounding corona. 

Future work will focus on exploring other corona 
geometries and applying the Monte Carlo scattering 
code to other QPO models in order to better under- 
stand the emission processes that describe the SPL 
spectral state. 
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